home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dpttrf.z / dpttrf
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPTTTTTTTTRRRRFFFF((((3333FFFF))))                                                          DDDDPPPPTTTTTTTTRRRRFFFF((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPTTRF - compute the factorization of a real symmetric positive definite
  10.      tridiagonal matrix A
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DPTTRF( N, D, E, INFO )
  14.  
  15.          INTEGER        INFO, N
  16.  
  17.          DOUBLE         PRECISION D( * ), E( * )
  18.  
  19. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  20.      DPTTRF computes the factorization of a real symmetric positive definite
  21.      tridiagonal matrix A.
  22.  
  23.      If the subdiagonal elements of A are supplied in the array E, the
  24.      factorization has the form A = L*D*L**T, where D is diagonal and L is
  25.      unit lower bidiagonal; if the superdiagonal elements of A are supplied,
  26.      it has the form A = U**T*D*U, where U is unit upper bidiagonal.  (The two
  27.      forms are equivalent if A is real.)
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      N       (input) INTEGER
  32.              The order of the matrix A.  N >= 0.
  33.  
  34.      D       (input/output) DOUBLE PRECISION array, dimension (N)
  35.              On entry, the n diagonal elements of the tridiagonal matrix A.
  36.              On exit, the n diagonal elements of the diagonal matrix D from
  37.              the L*D*L**T factorization of A.
  38.  
  39.      E       (input/output) DOUBLE PRECISION array, dimension (N-1)
  40.              On entry, the (n-1) off-diagonal elements of the tridiagonal
  41.              matrix A.  On exit, the (n-1) off-diagonal elements of the unit
  42.              bidiagonal factor L or U from the factorization of A.
  43.  
  44.      INFO    (output) INTEGER
  45.              = 0:  successful exit
  46.              < 0:  if INFO = -i, the i-th argument had an illegal value
  47.              > 0:  if INFO = i, the leading minor of order i is not positive
  48.              definite; if i < N, the factorization could not be completed,
  49.              while if i = N, the factorization was completed, but D(N) = 0.
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.